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Abstract 

The perturbation method in supersymmetric quantum mechanics (SUSYQM) is used 
to study the spheroidal wave functions' eigenvalue problem. Expanding the super-potential 
in series of the parameter a, the first order term of ground eigen-value and the eigen- 
function are gotten. In the paper, the very excellent results are that all the first two 
terms approximation on eigenfunctions obtained are in closed form. They give useful in- 
formation for the involved physical problems in application of spheroidal wave functions. 
PACs:11.30Pb; 04.25Nx; 04.70-s 



Introduction 



Since 1930s, spheroidal wave functions have made strong contributions to extensively 
theoretical and practical applications in pure mathematics, applied mathematics, physics 
and engineering. They appear in the fields, such as wavelet, random matrix, non-commute 
geometry, gravitational wave detection, quantum field theory in curved space-time, black 
hole stable problem; 3G mobile and broad band satellite telecommunication; steady flow 
of a viscous fluid and so on[I|-[S]. Nevertheless, they perhaps are one of the hardest work 
for researchers. They welcome the new thought and methods to deal with them. Since the 
appearance of the supersymmetric quantum mechanics (SUSYQM), its great power of solving 
the differential equation attracts tremendous attention. 

In this paper, we first use SUSYQM to study The spheroidal differential equations. First, 
brief introduction to the spheroidal problems and the ordinary methods to treat them. The 
spheroidal differential equations are 



dx 



(l-x2) 



d_ 

dx 



+ E + ax^ 



m 



x^ 



Q = o,xe (-i,+i). 



With the condition Q is finite at the boundaries x = ±1, they consist of the Sturm-Liouville 



eigenvalue problem. The parameter E can only takes the values Eq, Ei 



, En 



which are 
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called the eigenvalues of the problem, and the corresponding solutions (the eigenfunctions) 
Oo, Oi, . . . , &n, ■ ■ ■ are called the spheroidal wave functions [I]-[3]. 

The equations ([1]) have two parameters: m, a. When a = 0, the spheroidal wave functions 
reduce to the spherical wave functions, that is, the associated Legendre- functions Pl^{x). 
Though the spheroidal wave equations are extension of the ordinary spherical wave func- 
tions equations, the difference between this two kinds of wave equations are far greater than 
their similarity [T]. The spherical wave equations belong to the case of the confluent super- 
geometrical equations with one regular and one irregular singularities, whereas the spheroidal 
wave equations are the confluent Heun equations containing two regular and one irregular 
singularities. The extra singularity makes it extremely difficult to solve them[l]-[3], so very 
little information have been obtained concerning the analytically exact solutions. Therefore 
approximate and numerical methods are two main resources to rely on for the problem. 

Traditionally, the three-term recurrence relation methods are used to evaluate the eigen- 
values and eigenfunctions of spheroidal wave functions: one could solve the transcendental 
equation in continued fraction form or its equivalent or by power series expansion etc[T]- 
[3], [6]- [13]. For the details of these methods and their advantage and disadvantage, one could 
see the reference |13j . These methods mainly work for the numerical purpose, and also rely 
heavily on the numerical method. However, all previous works concentrate largely on the 
calculations of the eigenvalues, and particularly emphasize the small parameter approxima- 
tion and large parameter limits form of the eigenvalues. Little effort has been devoted to the 
related eigenfunctions due to the difficulty and complexity. 

Here, we give brief review on the eigenfunctions in small parameter approximation. No 
matter what method may be used, all eigenfunctions come into this kind of form in the end 

oo 

e„(x) = /'„"^ + 5]G„,(x)a«. (2) 

q=l 

Though there are many excellent works on the eigenvalues' approximation of small param- 
eter, no good works exist for the eigenfunctions approximation B„g(x). Even for the first 
eigenfunctions approximation, the existing works are only in the series' form 

e„i(x)= J2 BrP::+rix) (3) 

with infinite numbers needed to evaluate. Obviously, this series form does not reveal much 
information about the eigen- functions, even about its 1st order approximation eigenfunctions 
Qni{,x) itself. 

In recent years, super symmetric quantum mechanics have attracted tremendous attention 
for solvable potential problems. They not only provide clear insight into the factorization 
method of Infeld and Hull [16] , but also greatly improve the methods to solve the differential 
equations. See reference [13] for review on its development. 

The spheroidal eigenvalue problem is treated by the method of SUSYQM. This is the first 
time for researchers to use SUSYQM to study the ground eigenvalue and eigenfunction (that 
is Oq) of spheroidal wave functions in the small parameter a approximation. 

In usual small parameter approximation method, the key concept is the eigenfunctions 
and they are expanded in the form ([5]). On the contrary, the super-potential is the central 
concept in SUSYQM, and is expanded in the series form of the parameter a. This new method 
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is applied to study the spheroidal equations and unexpected results are obtained: the ground 
eigenfunction of the first order is in closed form, this in turn gives useful information on the 
eigenfunction and is helpful for their application. We also get the ground function for higher 
order terms in parameter a. 

The ground eigenvalue and eigenfunction in the first order 



In the following, we will use the new perturbation method in supersymmetry quantum to 
resolve the spheroidal eigenfunctions' problem. 

Though the form ([T]) is more familiar for research, the problem is easier to solve in the 
original differential equation than in the equation ([1]) . The original form is obtained from the 
eq.dH) by the transformation 

X = cos 0, (4) 



that is, 



1 d 



sm ( 



de 



+ a cos 



m 



sm 



e 



-Ee 



(5) 



the corresponding boundary conditions become G is finite at 9 = 0, vr. 

One writes the eqn.([5]) in the form of the Schrodinger equation by the transformation 



G 



the differential equations turn out to be 



+ 



1 2 

- + a cos 
4 



sin 2 



m 



sm 



+ E 



^' = 



and the boundary conditions become 

From the equation ([7]), one knows the potential is 



(6) 



(7) 



(8) 



V{0, a, m) 



1 2 

a cos 

4 



+ 



m 



sm 



(9) 



The super-potential W is determined by the potential V{0^ a, m) through the Reccita's equa- 
tion 



W' = V{0, a) - Eo 



(10) 



where the substraction of the ground energy just makes the eqn.d?]) factorable. Actually this 
equation is the same hard to treat as that in the original form ([7]). The approximate method 
naturally comes to one's mind. Hence, when the absolute value of a is small, it is the the 
super-potential W that could be expanded as series of the parameter a, that is. 



W = Wo + aWi + a^W2 + a^Wa + .... 



(11) 
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_ ^/ ^ ^,2 _wl^ + a {2WoWi - Wi) + (2W0W2 + Wf - W^) 
+ {2W0W3 + 2W1W2 - W^) + {2W0W4 + 21^1^3 + TV| - Wi) + . . . . (12) 
One can write the perturbation equation as 



°° 1 ^2 _ 1 00 

-W = V{e, a, m) - y 2Eona" = -t + - a cos^ - V 2Sona" (13) 

n 4 sin 9 ^ 

ra=0 n=0 

There are two lower indices in the parameter E^n with the index refereing to the ground 
state and the other index n meaning the nth term in parameter a. The last term J2'n=o ^-Eono" 
is subtracted from the above equation in order to make the ground state energy actually zero 
for the application of the theory of SUSYQM. Later, one must add the term to our calculated 
eigen-energy. Comparing the equations (fT2]) . Q, and (fT3|l . one could get 

1 77^2 _ i_ 

Wi -Wl, = -- + - 2Eoo (14) 
4 sm 

2W0W1 - W[ = -a cos^ e - 2Eqi (15) 

2WqW2 + Wl -W2 = -2Eo2 (16) 

21^0^3 + 2W1W2 - VF3 = -2^03 (17) 

2W0W4 + 21^1^3 + Wi -Wi = -2^04 (18) 



(19) 



From the eq. p^ . we get 



Wo = -(^m + ^^ cot 9, 2Eoo = m{m + 1). (20) 

Then, we can write the other equations more concisely 

P^i + (2m + 1) cot 6*1^1 = cos^ 6* + 2^01 (21) 

W2 + (2m + 1) cot eW2 = Wf + 2Eq2 (22) 

W3 + (2m + 1) cot 6*^3 = 2W"iW2 + 2-^03 (23) 

W4 + (2m + 1) cot = 2WiWz + + 2^04 (24) 

W5 + (2m + 1) cot eW^ = 2WiWa + 2W2W^3 + 2^05 (25) 

; (26) 

After obtaining the zero term Wq for the super-potential W , the first order Wi can be gotten 
as 



with 

"^^(cos^ + 2^01) (28) 



'^^^ - sin^--^ 



d9 



Ai = j i^W"" 9 cos^ 9 + 2Eoism^'^ 9) 9 dO (29) 
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Suitably changing the independent variable to x = cos and expanding the term as following 



sin^-^ = (l-cos'ey =(l-x^Y = Y,{-lf 

where (X) = ^x(m-i)x^..x(m-fc+i) ^ ^^^^ ^^^^^^^ 

m 

So, 



m 



„2k 



k=0 



^1 



m ' 
, k I 



k=0 



008"^^+^ 9 2^01 cos^'^+i 
+ 



2k + 3 



2k + 1 



Er=o(-i)'+nT) 



2k+3 



2k+l 



(30) 



(31) 



(32) 



The quantity £?oi needs to be determined by the boundary conditions. This in turn require 
to calculate the ground eigenfunction upon to the first order by 



A^exp 
iVexp 



wde 



N sin"^+ 2 exp 



J Wod9 - a J Wide 
a J Wide 



* expO(Q;^) 
* expO(a^). 



(33) 
(34) 

(35) 



Whenever the eigenfunction is obtained, the boundary conditions ^'|6i=o = ^|e=7r = would 
choose the proper Eqi. This sounds very easy, but is a tough task in reality. The complete 
calculating process is left to the appendix 1. The results are 

2Eoi = -^^-^ (36) 



Wi 



2m + 3 
sin e cos 9 



(37) 



2m + 3 

With the first order term of the super-potential Wi, one could compute the second term 
W2 by the same process. This is not easy either. The appendix 2 gives the detail of the 
calculation. Nevertheless, the results arc elegant: 

2m + 2 



2Eo2 = 



W2 



(2m + 3)3(2m + 5) 
— sin cos 9 

+ 



sm f/ cos ( 



(2m + 3)3(2m + 5) (2m + 3)2(2m + 5) 
The ground eigenfunction upon to the second order becomes 



N exp 
N exp 



Wd9 



- J Wod9 - a J Wid9 - j W2d9 



expO(a^) 



(38) 
(39) 

(40) 

(41) 



= (sine)'^^^^ exp 



a sm u 
4m + 6 



* exp 



sin^ ( 



sin^ I 



2(2m + 3)3(2m + 5) 4(2m + 3)2(2m + 5) 



* expO(Q;^). 



(42) 
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When m = 0, These results p6|) . (|37|) . (j38|) . (|39|) reduce respectively to 

2Eoi = -L Wi = ^ sin 6* cos 6*; (43) 

2^02 = -t|-, H/'a = -t^ sin cos + ^ sin^ ^ cos 9. (44) 
loo loo 45 



They are in complete accordance with that in reference [18] where the spheroidal wave 
functions are treated by SUSYQM in the case m = 0. 

Back to equation ([T|) with their relationship ([6]), we could obtain 

a(l- 



2 

exp 




2\2 



a'{l-x') a'{l-x') 



""""P 2(2m + 3)3(2m + 5) 4(2m + 3)2(2m + 5) *«^P0("')- (^5) 
Expanding the exponential functions in the above equation, the results elegantly turn out 

as 

^2 / 1 ^2 n ^2\2 



2/ 1-^' (1-X^)^ 

1-"-; ;: + a hrTT^ ^t^ttt^ ^ + 



4m + 6 \^2(2m + 3)3(2m + 5) 8(2m + 3)(2m + 5)^ 

+0(a=^). (46) 

The above equation clearly shows how the ground state function changes as the function 
of variable x when a is small and could be compared to the non-perturbation case = 
(1 — x^)~. The result is much better than the usual result of the series form 

oo oo 

Qoix) =P^ + aY, i?,C^+,(x) + a'Y. + 0{a^) (47) 

q=l q=l 

with infinite numbers Bq, q = 1, 2, ... need to be determined for the first order term and 
Cq, q = 1, 2, ... need to be determined for the second order term. Prom another point 
view, the results in the eqns. (|45| ). (j46|) give the method to determine the infinite numbers 
Bq, Cq, q = 1, 2, .... 

In conclusion, SUSYQM provides a new opportunity to treat the spheroidal wave functions 
and indeed they give new results in the eqns. ()45p . ()46p . Further calculations can be done by 
the same way, nevertheless, the higher order term Wn is more complex than the lower one. 
The maximus of Wi and W2 satisfy 



1 y 11 1 

max 1^1 = -, maxW2 = ttt^ < < maxW'i. (48) 
D 21dU 54U 

Though the further calculation is not processed here, the reasonable guess is that the higher 

term Wn is, the smaller its maximus. The guess mainly comes from the calculation of the 

quantity Wi, W2, see the appendix 1 and appendix 2 for details. If the guess is right, 

00 

wie) = J2 (49) 

n=0 

might be analytic function in all complex plane a, the only singularity of the function W 
as the variable a is a = 00. Only the calculation in the appendixes can process on and on, 
could the guess be tested. Might some day the computer can do the work, this is the reason 
that the calculations in appendixes are extremely detailed. 
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Appendixl: Simplification of the super potential of the first 
term Wi 



(50) 



The calculations are very complex, we rewrite the eq. ()32p here again for convenience 

Er=o(-i)'+HT) 



2k+3 



2k+l 



sm 



2m+l , 



As state before, the quantity £"01 is determined by the requirement that the eigenfunction 
is zero at the boundaries 9 = 0, tt; this in turn demands / Wid9 finite at the boundary. 
Therefore, the calculation of the term J Wid9 is first processed . By transformation 



r = sm f 



(51) 



and denoting / WidO by / , it reads 

Er=o(-i)'+HT) 



Wid9 



2fc+3 2fc+l 



p2m+l 



dr 



Using formula 



, , 1 k+1 / 

(i-rf« = i:(-iy( 

1=0 \ 



(52) 
(53) 



it becomes 

/m 

k=0 

m 
k=0 
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m 
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k 





Et^'(-i)'rTV 



k+l\9l-2m-l 



2k + 3 



+ 



2Eo,Ef=oi-^y{i )r 
2k + 1 



l(k\^2l-2m-l 



(_l)/(fc+y2^^2m k 2£;o,(-i)^(fc)r^^-2m 
^ (2A; + 3)(2/ - 2m) (2A; + 1)(2/ - 2m) 



EE 

k=0 1=0 



(2fc + 3)(2/ - 2m) 



+ EE 

k=0 1=0 



{2k + 1)(2; - 2m) 



dT 
(54) 
(55) 



In order to exchange the sums order in the above equations, one must notice the fact that 

and I < k ^ k > I 



, , -, (^\k + l + l(m\(k + l\ 21-27, 

l<k + l ^ /fc > / - 1 in the first term ^^=0 Ef=o — ' 



(2fc+3)(2Z-2m) 
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in the second term ^^=0 Eto ^ (i+ijiiAj 



So 



m+l m 

E E ■ 

Z=0 fe=/-l>0 



-1) 



I 



k+l+l(m^^(k+l\2l-2m m m 

+ EE 

1=0 k=l 



2i?oi(-l)*^+'+nT)(z)^""'"' 



E 

1=0 



(2A; + 3)(2/ -2m) 



(2/c + 1)(2/ - 2m) 



,,=^i>o + 3)(2/ - 2m) + ^ {2k + 1){21 - 2m) 



T + 



(56) 
4m + 6 



where the first sum under the condition I = m + 1 becomes 



Defining 



it gets 



4m+6 ■ 



iVi 



2 I 



-1>0 



(2A; + 3)(2/ - 2m) 



Z = 0,l,2, 



V ^ ^> 7 = 012 

^^(2fc + l)(2/-2m)' ^ 



(58) 
(59) 



J^[N,i + 2E01N2 i] r2'-2- + 



1=0 



4m + 6 



(60) 



There are terms the above equation. With the eq. ()5ip and the fact I < m, 

these term as sin^'~^™ become infinite as ^ ^ 0, vr. By eq. (|35p and the eigenfunction's 
boundary condition at = 0, vr, the coefficients of those terms must be zero. There is only 
one quantity £"01 unfixed, could one choose proper £"01 to to make the eigenfunction finite at 
the boundaries? Actually, one only has one choice to select £01 by 



iVi + 2£oiA^2 = 0, Ni 



j_vLir!(T) ;v 

2m (2A: + 3) ' '''' 



1 ™ ( — ^'\k+2(m\ 

E \oi , (61) 



2m 



fc=o 



{2k + 1) 



do the other terms in eqn. (j60|) automatically become zero under the condition (j61l ). 

Fortunately, this can be done and the following is the proof. The inductive reasoning is 
used to give the proof. In order to determine the quantity £01 under the condition (|6ip . one 
must first simplify A'^i O) -^2 0- From the formula ( see reference [ISl on page 389) 



(1 - T^rdT 



' cos^^+i ed9 



(2m) 



and the similarly one 



f\\i-T^rdT-- 

Jo 

(2m)!! (2m + 2)!! 



1 - r'^TdT 



(2m + 1)!!' 



:i - r'^T+^dT 



(2m + 1)!! (2m + 3)!! 



(62) 

(63) 
(64) 
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and also with the formula (|53|) . it is easy to obtain 



(1 - r^rdr = / E T (-1)'-"^- = E ^r-K (66) 



t n \ ^ / t n 2A: + 1 

fe=0 \ / fe=0 



Comparing these equations with that of ()6ip . one reaches 

AT = ^v(-l)'^'(T) _ 1 r (2m)!! (2m + 2)!! 



2m 2k + 3 2m 

k=0 



(2m + 1)!! (2m + 3)!! 



AT = 1 f. (-l)^+HT) _ 1 (M!! .g7^ 
^° 2m 2/c + l 2m(2m + l)!!' ^ ^ 

Consequently, the quantity 2Eqi is simplified as 

2Eo^ = -^ = —. (68) 

A^2 2m + 3 ^ ^ 

A''i + 2E01N2 = is guaranteed by the choice in eqn. ([68]) . According to inductive 
reasoning, one needs to prove 

Ml ;+i = iVi i+i + 2E01N2 i+i = (69) 

under the assumption that 

Mii = Nii + 2E01N2 I = 0. (70) 

The key idea is to find the connection between the terms Mi i = Ni i + 2£^oi-^2 i and 
Ml ;+i = A^i i+i + 2E01N2 i+i. By the definitions of iVi i, N2 i, I = 0, I, 2, . . . in eqns.(l58])- 
([59|) . one has 

"1 (_-\\k+l+2(m\(k+l\ 

- ^(2A: + 3)(2Z + 2-2m)' 

AT - V ^) \k)\l+l) /^^.N 

^^'^^ - ,E^(2fc + l)(2/ + 2-2m)- 



Due to the following relation 



A + ^ (fc + 1)! ^ (fc + 1)! 



_ k+l-l 

2k + 3\l + lJ {2k + 3){k-l)l{l + iy. {2k + 3){k + l-l)W.* (l + 1) 

(fc + 1)! ^ fe + |-(/ + l) 
(fc + 1 -/)!/! * (Z + l)(2A; + 3) 

/ + 1 2k + 3\ I J 2l + 2\ I y' 



(73) 
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it is easy to get 

{21 + 2- 2m)Ni i+i 



(-1)^+'(T)(S) 



E 

k=l 
m 

E(-i)' 



2fc + 3 



(74) 



k+l 



k=l 



m ' 
, k , 



( 



k + V 



+ 



1 (k + V 



l + h 



l + l 2k + 3\ I J 2l + 2\ I 



(75) 



= E 

k=l-l>0 L 
m 

+ E 



1 



l + l 2k + 3 
(-1)^+' (m\(k + l 



■1) 



rn \ (k + l 



2(1 + 1) 



k=l 



21 + 2 \k 



Now, the calculation becomes 



{21 + 2- 2m)N 



1 l+l 



{21 - 2m)Ni I 



I 



(76) 



m 



{I + 1) , 
under with the help of eqn. (|58p . Definition of the quantity 



1 (k + l^ 



Qii 



ii-d , g^ (-i)^+'(T)('r) 



21 + 2 



k=l 



21 + 2 



I 



(78) 



may make the equation (|77p simple as 



{21 + 2- 2m)Ni 



l + \ 
G + 



{21 - 2m)Ni i + Qiil = 0, 1, 2, . . . , m-1. 



(79) 



In completely similar way, one could get the relation between A^2, i and N2 i+i- the 
formula 



1 



k 



1 



k\ 



1 



kl 



k-l 



2k + l\l + l 



2k + l{k-l- 1)\{1 + 1)\ 2k + l{k- /)!(/)! 

k\ {k + \)-{l + \) 

{k-l)\{l)\* (2A; + !)(/ + !) 

l + h 



1 



l + l 2k + l\l 
helps to simplify the quantity N2 i+i 



k\ 1 k' 

+ 



2l + 2\l 



(80) 



{21 + 2- 2m)N2 l+l 



E 

k=l+l 



2k + 1 



E (-1) 

k=l+l 



k+l 



m 



1 + ^ 1 ( k\ 1 ( k 
l + l ^ 2k + l\l ) 2/ + 2 1^/ 



(81) 
(82) 
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f^ll + l^ 2k + I \k)\l) l + l^ 21 + 1 
Using the eqn. (f59]) and the fact — x '^-^^ — (T)(D ~ 2lT2(T)' ^^^^ *° obtain 
The similar definition of the quantity 



(85) 



makes the equation ()84p simply become 



(/ + i)(2/-2m) , , 

(2/ + 2 - 2m)iV2 /+i = ^ « + Q2 ^ = 0, 1, 2, . . . , m - 1. (86) 

The eqns.(I79D,(l86D teh that 

Mi;+i = A^i i+i + 2£oiiV2 /+/ (87) 
(/ + i)(2/-2m) 

= (i Vl)(2i -2,. + 2) ^^' . + [Cm + 2E„Q, ,1 . I = 0, 1, 2, . . . . ™ - 1. 

the relation between Mi ;,-/Vfi would be the desired result if one could prove that Q\ i 
Q2 I = 0, 1 = 0, 1, 2, . . . , m — 1. This is the following work. 
Using the formula 

m / \ 

under the special case of / = r, g = {1 — r)™, one gets 

[ril-rrf^ = T[{l-Trf^ + (^^[{l-rrf-'^'^ (90) 

= (-l)'m(m-l)...(m-/ + l)T(l-r)"-' 

+ (-l)'-i/ xm(m-l)...(m-/ + 2)(l-r)™-'+i. (91) 



The calculation may also go by diff'erent way: 



[r(l - rn 



E(-i)' 

.k=0 



m 

. k , 



(0 



(92) 
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m 



k=l-i 
1 



{k + l)k{k - 1) . . . {k - I + 2)t''-^+^ (93) 
k I \ (k + \ 



k=l-l 



k 



I 



(94) 



Comparing the results of eqns. ()9ip . ()94p and taking the function [r(l — r)™"]^'^ at r = 1 
under the condition / < m, the good result reaches 



E (-1) 

fc=Z-l 



m\ fk + 1 



k 



I 



T 



k-l+l 



I 



T=l 



By eqn. ([78|) . it is easy to get 

Qii = 0, / = 0, 1, 2, 
Similarly, one also could prove that 

Q21 = 0, 1 = 0, I, 2, 

Therefore, the eqns. 



under the condition 



m — 1. 



m — 1. 



0. 



(I96I1 . (1971) imply that 
Ml z+i = 0, / = 0, 1, 2, 



m — 1 



by induction one gets 



Mi = 0; 



M„ = 0, n = 1,2,... ,m. 



(95) 



(96) 



(97) 



(98) 



(99) 



(100) 



Hence, the boundary conditions could be satisfied by just selecting the only one quantity 
2^01 = -2^- With the good results Mn = N11 + 2E01N2 i = 0,l = 0, 1, 2, . . . , m , one 
can greatly simplify the first order super-potential Wi in the eqn. (j50|) . By similar method as 
before, rewrite Wi by changing the independent variable to 



r = sm f 



(101) 



and expanding terms (1 — t^)'^+^, (1 — r^)'^, that is 



Er=o(-i) 



Wi = 

m 

(l-r^)^E(-l 

k=0 



2k+3 



2k+l 



T 



2m+l 



Ef=V(-i)HT). + i^oiEf=o(-i 



HI) 



2k + 3 



2k + 1 



(102) 
'2'-2"^-\l03) 
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then changing the order of the sums just as before and dividing Wi by (1 — r^)2 for good 
looking in the formula, it reads 



(l-r2)2 



m+1 rn ^_'^^k+l+l ^k+1^^2l~2m-l 

/=0 fc=«-l>0 



+ EE 

/=0 fc=/ 



(2A; + 3) 
(2A; + 1) 



(104) 



m m 



( _-\\k+l+l (■rri\ 
^ ^ V -^^ Vfc/W I ^2l-2m-l 



+ 



1=0 k=l-l>0 



{2k + 3) 



.2i-2m-l 



(2fc + 3) 



/ -i\k+l+l(m\(k\ 
\ ^) \k)u) ^21-2: 



l=m+l 
k=m 



E 

/=0 



V (-i)^+'+HT)('r) , op 

E PT3) + 2^oiE 



fc=i-l>0 



{2k + 1) 



^2l-2m-\ 



+ 



2m + 3 

[iVi i + 2^01 i] * {21 - 2m)r2'-2"^-i 
r 



+ 



2m + 3 



2m + 3 

So the quantity Wi could be written tidily as 

r(l-r2)l 



2m + 3 



sm a cos b 
2m + 3 



(106) 
(107) 
(108) 

(109) 



Appendix2: Simplification of the super potential of the first 
term W2 

Though the process of the calculation of W2 repeats that of Wi in the appendix 1, there 
still are some needs to write it down. Perhaps this may be used for further computation, 
even by computer. 



where 
A2 



J (W^f + 2£:o2)sin 

/ 



W2 



ede 



sm 



2m+l , 



(2m + 3)2 
1 



gij^2m+3 Q ^Qg2 Q ^ 2^02 Sin2"*+1 9 



(2m + 3)2 



{l-x^r+^x'' + 2Eo2{l-x 



2\m 



de 

dx 



. . . ,X = COS ( 



(110) 

(111) 
(112) 
(113) 
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.2k 



m+1 



-1) 



k+lfm+l\ 2k+3 



(2m + 3)2(2fc 

m+l ^_'^'^k+l fm+l\ 
k=0 



(2m + 3)2(2A; + 3) 
The other important term is 

11 = J W2de 



+ 2Eo2 



fe=0 



2k + I 



m+l 



^ (2m + 3)2 (2A; + 3) sin^^^+i °^ ^ (2A; + 1) sin^^^+i 9 



dx{lU) 
(115) 
(116) 

(117) 
de. (118) 



By the transformation r = sin0, one may have 

j Wide = 



y (-l)fe+i(-+i)(l-r2)^'+i - (_i)^+i(™)(i_^2)fc 

^ (2m + 3)2(2fe + 3)T2'"+i 



fc=0 



fc=0 



(2/c + l)r 



2m+l 



m+l fc+l /_i xfe+Z+l/m+lw/c+lN 2«~2m m k ( 



A;+i+l('m\ (k\ 2l-2m 



-Oil) 



k=0 1=0 



(2m + 3)2(2/e + 3)(2/-2m) 



A;=0 i=0 



(2A; + l)(2l-2m) 



<il9) 



(120) 



Exchanging the sum order, it reads 



m+2 fc+1 



/_1\fc+i+l('m+l\/fc+l\ 2i-2m m m 

^-^ UiUil + 2£;o2EE 



(__l)fe+i+l(™) l^^j.2l-2m 



(2m + 3)2(2A; + 3)(2Z-2m) 



i=0 fe=i-l>0 

^ ^ (2m + 3)2(2/c + 3)(2/-2m) 



/=0 k=l 

rn rn 



{2k + 1){21 - 2m) 



_ 1 \k+l+l (m\ (k\ 2l-2m 

+ 2^02 E E ' L , .wo, (122) 



i=0 fc=i-l>0 
m+2 m+l 

^.=?+i.=S>o (2"^ + 3)^(2fc + mi - 2m) 



Z=0 k=l 



{2k + 1){21 - 2m) 



E 



m+l ^_-|^-jfe+/+l ^m+lj ^fc+l-j^2«-2m m -[^-jfc+i+1 ^mj ^fcj^2i-2m 

^^^^^ (2m + 3)2 (2fc + 3) (2/ - 2m) + 



(2fc + l)(2i - 2m) 



m+2 m+l ^_-|^-jfc+«+l ^m+l-j ^fc+l-j^2«-2m 

^ .=?+! .=S>o (2m + 3)2(2fe + 3)(2/-2m) 
It is better to define the following quantities 

(_l).+/+l(m+l)(.+l) 



(123) 
(124) 
(125) 



AT, 



3, / 



m+l 

E 

k=l-l>0 



(2m + 3)2(2fc + 3)(2l - 2m)' 



/ = 0, 1, 2, 



m 



(126) 



™ ^_^\k+l+l(m\(k\ 



k=l 



{2k + 1){21 - 2m) 



, 1 = 0, 1, 2, . . . , m 



(127) 
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and simplify the last term 



m+2 m+1 

E E 



i=m+l fe=/-l>0 

(2m + 6)r 



(2m + 3)2(2/? + 3) (2/ -2m) 

^3 



+ 



(2m + 3)3(2m + 5) 3(2m + 3)2(2m + 7) ' 



Then the express becomes 



II = J2[N3l + 2^02 iV4 l] T 



2/-2m 



(2m + 6)T 



+ 



1=0 



(2m + 3)3(2m + 5) 3(2m + 3)2(2m + 7) 



(128) 
(129) 

(130) 



Just as done before, one needs to select proper £"02 to make all / + 2£'o2A''4 i zero. By com- 
paring the quantities A'^i i,N2 uN-^ i, N4 i, the good relations among them could be revealed: 



N4i = N2i, 1 = 0, I, 2, m 



(131) 



and the great similarity between A'^i i,N-3 i. The transformation forms between A'^i i-^-i,Ni i, 
l^ijN^ I are similar. 



3, l+l 



m+1 

E 

k=l>0 



(2m + 3)2(2/c + 3)(2/ - 2m + 2) 



By the relation ([73 



2k + 3\l + lJ~ l + l^ 2k + 3\ I J^2l + 2\ I 



one could get 



m+1 (_-]\k+l(m+l\/k+l\ 

(2m + 3)2(2/ + 2 - 2m)iV3 = ^ ^ ^ ^ k 

k=l 



2k + 3 



m+1 

E(-i) 

k=l 
m+1 

E 

k=l-l>0 
m+1 



k+l 



m+1 

, k , 



1 + 



1 



1 A + 1 /A; + l' 

/ + 1 ^ 2fc + 3 \ / J ^ 2/ + 2 i / , 



/ + i 



1 



/ + 1 2A; + 3 



(-1) 



k+l- 



_^/m+l\ (k+l 



/m+l\ 

V ;-i / 
2(/ + l) 



+ E 



(-1) 



k+l 



m + l\ k+l 



k=l 



21 + 2 



with the help of eqn. p26|) . the calculation becomes 

(T-Y) 



(;+i)(2/-2m) 
- (/ + l)(2/ + 2-2m)^=^ ' + 



/-m + lWfc+lN 

V i-1 / _)_ Y^m+1 / 1 \ fc+i I fc A i J 



2i+2 



(2m + 3)2(2/ + 2-2m) 



(132) 



(133) 



(134) 
(135) 



(136) 



(137) 
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The definition of the quantity 
Q31 = 



21 + 2 



+ E 

k=l 



21 + 2 



21 + 2 

makes the above equation ()137p become 
{l + \){2l-2m) 



+2.5^1 ^ V ^ ; 



m+l\ fk + l 



(/ + l)(2/ + 2-2m)' 
Another useful relation is 



N31 + 



Q31 



4 i+i 



(2m + 3)2(2/ + 2-2m)' 
(^ + i)(2/-2m) 



1 = 0, 1, 2, 



2 i+i 



(/ + l)(2/ + 2-2m) 



2 z 



(/ + i)(2/ - 2m) 



4 / 



(/ + l)(2/ + 2-2m) 
by the use of Q2 i = 0. Now the quantities M2 i is defined as 

M2i = N3i + 2EQ2Niu / = 0, 1, 2 
then the relation between M2 i and M2 i+i is 



, . . . , 



ra 



(Z + i)(2/-2m) 



2 «+i 



Mo ,+ 



3 I 



{I + 1)(2/ - 2m + 2) ^'^^ ' ' (2m + 3)2(2/ + 2 - 2m) 
In the following, it will be proven that Q3 i = 0- With the help 



k=0 



under the special case f = t, g = {1 — t)"^~^ , one gets 



t(1 - T 



\m+l 



(0 



(1-r) 



m+l 



(0 



+ 



(1-r) 



m+l 



(l-l) 



(138) 

m- 1.(139) 

(140) 
(141) 

(142) 



/ = 0, l,...,m-l. (143) 



= (-l)'m(m + 1) . . . (m - / + 2)r(l - t)™-'+^ 

+ (-1)'"^/ X m(m + 1) . . . (m - / + 2)(1 - r)™-'+2. 

The alternative way to compute is 



t(1 - T 



(0 



m+l 

E(-i) 

Lfc=0 
m+l 

E(-i) 

fc=Z-l 



/ + 1 

■ k 



1^1 m + l 



n(0 



^fe+1 



(144) 

(145) 
(146) 

(147) 



{k + l)k{k -l)...{k-l + 2)t'=-'+i (148) 

(149) 



^^m+1^ ( k + 1 
I 
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The valuation of [r(l — 7-)™+^] '•'^ at r = 1 under the condition I < m + 1 



IS 



m+l 

k=i-i 



m + l\ k + l 



^k-l+i 



r=l ^=='-1 



From the eqn. ()138p . it is easy to get 

Q3« = 0/ = 0, 1, 2, m-1. (151) 
Finally, one gets the good result 

So the choice of M2 = guarantees M2, ; = 1, 2, . . . , m. The quantity £'02 is obtained by 

Eo2 = -^. (153) 

iV4 



As in the case A'"! 



^30 = 



m+l ^_l^k+2^m+l^ 



2m ^ (2m + 3)2(2A; + 3) 2m(2m + 3) 



(2m + 2)!! (2m + 4)!! 
(2m + 3)!! ~ (2m + 5)!! 



- ^^°-2;^^^^ 2/fc + l -2;^(2m + l)!!- ^^^^^ 



Hence 



The expression 



_ A^s _ 2m + 2 

°' " " A4^ " "(2m + 3)3(2m + 5)- ^^^^^ 



Z^fe=0 (2m+3)^(2fc+3) ^ ^-^^02 Z^fc=0 2fe+l 
^2 = . 2m+l/) (156) 

sm ^ y 

could greatly simplified by r = sin0 and the use of the elegant formula M2 ; = 0: 

W2 y (-i)^+^(-+^)(i-t2)^+i f. (-i)^+HT)a-^T 

cose (2m + 3)2(2A: + 3)r2"^+i °^ (2A: + 1)t2™+i ^ ' 

m+l fc+l \A:+i+l ('m+lwfc+lN m k ( _^\k+l+l (m\ (k\ 

^ (2m + 3)2(2A; + 3)r2"^-2'+i ^ A. (2A; + l)r2"^-2/+i ' 
Exchanging the sums order, it is easy to get 

m+2 m+l / 1 Nfc+i+l/m+lN /fc+l\ m m / n fc+i+i /m\ /fe\ 

COS e ^_fl^n (2"^ + 3)2 (2A; + 3)r2™-2Z+i ^ ^^02 2^ (2A; + l)r2™-2«+i^''^^^ 
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comparing the eqns. (jl26p and (jl27p . one may have 



W2 YT=o [^3 / + 2^02 iV4 /: 



COS ( 



m+2 m+1 

+ E E 



(_l).+m(n.+l)(.|l) 



r2-2/+i ' ^J^^^ ^J-^^_^ (2m + 3)2 (2fc + 3)r 

(_l).+.+l(m+l)(/c+l^ 



m+2 m+1 

z=?+i.=S>o(2- + 3P(2^ + 3)-^'"-^'+^ 



+ 



(2m + 3)3(2m + 5) (2m + 3)2(2m + 5) 
by the use of M2 ^ = iVa ; + 2£'o2iV4 ; = 0. The elegant form of W2 is 



W^2 



-1 



+ 



sm 



(2m + 3)3(2m + 5) (2m + 3)2(2m + 5) 
This ends the appendix 2. 



sm u cos ( 



2m-2/+T (160) 

(161) 
(162) 

(163) 



References 

[1] C. Flammer 1956 Spheroidal wave functions. ( Stanford, CA: Stanford University Press). 

[2] J. Stratton et al., Spheroidal Wave Functions Wiley, New York, 1956 

[3] Le-wei Li, Xiao-kang Kang, Mook-seng Leong, Spheroidal wave functions in electromag- 
netic theory. (John Wiley and Sons, Inc., New York, 2002) 

[4] S. A. Teukolsky Phys. Rev. Lett. 29 1114 1972; S. A. Teukolsky 1 J. Astrophys. 185 635 
1973. W.H. Press and S. A. Teukolsky 1973 J. Astrophys. 185 649; B. White 1989 J. 
Math. Phys. 30 1301; J. M. Stewart 1975 Proc. R. Soc. Lond A 3U 65; J. B. Hartle and 
D. C. Wilkins 1974 38 47; R. A. Breuer, M. P. Ryan Jr and S.Waller 1977 Proc. R. Soc. 
Lond A 358 71; S. Chandrasekhar 1983 The mathematical theroy of black hole; (Oxford: 
Oxford University Press) ch8-ch9; E. W. Leaver 1986 J. Math. Phys. 27 1238. 

[5] D. Slepian and H. O. Pollak 1961 Bell. Syst. Tech. J 40 43; D. Slepian 1964 Bell. Syst. 
Tech. J 43 3009; H. Xiao, V Rokhlin and N Yarvin 2001 Inverse Problem, i7 805. 

[6] J. Caldwell, J. Phys. A 21, 3685 1988 

[7] E.T. Whittakar and G. N. Watson 1963 A Course of Modern Analysis (Cambridge 
University Press) chXVII p366. 

[8] D. B. Hodge, J. Math. Phys. 11, 2308 1970 

[9] B. P. Sinha and R. H. MacPhie, J. Math. Phys. 16, 2378 1975. 

[10] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, 9th ed. Dover, 
New York, 1972. 

[11] P. E. Falloon, P. C. Abbott, and J. B. Wang, math-ph/ 0212051. 



18 



[12] L.-W. Li, M.-S. Leong, T.-S. Yeo, P.-S. Kooi, and K.-Y. Tan, Phys. Rev. E 58, 6792 
1998. 

[13] E. Berti, V. Cardoso, M. Casals, Phys. Rev. D 73, 024013 (2006); E. Berti, V. Cardoso, 
K. D. Kokkotas, and H. Onozawa, Phys. Rev. D 68, 124018 (2003); E. Berti, V. Cardoso, 
and S. Yoshida, Phys. Rev. D 69, 124018 (2004); [20] E. Berti, g r-qc/041102"5[ 

[14] F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251, (1995) 268, and references 
therein. 

[15] R. Dutt, A. Khare, and U. Sukhatme, Am. J. Phys. 56, (1988) 163; A. Khare and U. 
Sukhatme, Jour. Phys. A21, (1988) L501; R. Dutt, A. Khare and U. Sukhatme, Phys. 
Lett. 181B, (1986) 295. 

[16] L. Infeld and T.E. Hull, Rev. Mod. Phys. 23 (1951) 21. 

[17] Guihua Tian, Shuquan Zhong, Arxiv: 0906.4685: Solve spheroidal wave functions by 
SUSY method, preprint 

[18] Guihua Tian, Shuquan Zhong, The ground state eigenfunction of the spheroidal wave 
functions, preprint 

[19] I.S. Gradsbteyn, L.M. Ryzbik, Table of integrals, series, and products. 6th edition, El- 
sevier(Singapore)pte. Ltd, 2000. 



19 



